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Abstract: Spin observables may reveal much deeper properties of non perturbative 

hadronic physics than unpolarized quantities. We discuss the polarization of hadrons 
produced in e + e~ annihilation at LEP. We show how final state qq interactions 
may give origin to non zero values of the off-diagonal element p 1 __ 1 of the helicity 
density matrix of vector mesons: some predictions are given for K*,<f),D* and B* 
in agreement with recent OPAL data. We also discuss the relative amount of vector 
and pseudovector meson states and the probability of helicity zero vector states. 
Similar measurements in other processes are suggested. 
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p x _ t (V) in the process e e+ — > qq — > V + X 

The spin properties of hadrons inclusively produced in high energy interactions 
are related to the fundamental properties of quarks and gluons and to their elemen- 
tary interactions in a much more subtle way than unpolarized quantities. They test 
unusual basic dynamical properties and reveal how the usual hadronization mod- 
els - successful in predicting unpolarized cross-sections - may not be adequate to 
describe spin effects, say the fragmentation of a polarized quark. 

We consider here the spin properties of hadrons created at LEP. It was pointed 
out in Refs. and that final state interactions between the q and q produced in 
e + e~ annihilations - usually neglected, but indeed necessary - might give origin to 
non zero spin observables which would otherwise be forced to vanish: the off-diagonal 
element p 1 _ l of the helicity density matrix of vector mesons may be sizeably different 
from zero |]J due to a coherent fragmentation process which takes into account qq 
interactions. The incoherent fragmentation of a single independent quark leads 
instead to zero values for such off-diagonal elements. 

We present predictions for p 1 _ 1 of several vector mesons V provided they are 
produced in two jet events, carry a large momentum or energy fraction z = 2E V / 'y/s, 
and have a small transverse momentum p T inside the jet. Our estimates are in 
agreement with the existing data and are crucially related both to the presence 
of final state interactions and to the Standard Model couplings of the elementary 
e~e + — > qq interaction. 

The helicity density matrix of a hadron h inclusively produced in the two jet 
event e~e + — > qq — > h + X can be written as [§], H 

1Sh g,X,X x ,X q ,X q ,X' q ,X' q 

where p(qq) is the helicity density matrix of the qq state created in the annihilation 
of the unpolarized e + and e~, 

P\,X^ q iQQ) = E M A A ;A_A + M KX-X_X + ■ ( 2 ) 

The M's are the helicity amplitudes for the e~e + — > qq process and the D's are the 
fragmentation amplitudes, i.e. the helicity amplitudes for the process qq — > h + X] 
the J2x,x x stands for the phase space integration and the sum over spins of all the 
unobserved particles, grouped into a state X. The normalization factors N h and 
N q q are given by: 

N h = E D \^\^ p\,\-^ (qq) d \^k^ = E D i > ( 3 ) 

q,X;X h ,X x ,X q ,X (j ,X' q ,X' lj g 
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where D% is the usual fragmentation function of quark q into hadron h, and 



*«=\ E I^V^aJ 2 - (4) 



The helicity density matrix for the qq state can be computed in the Standard 
Model and its non zero elements are given by 



l (g v - g 



a) 2 

p+-,+-(q9) = - 2 p + ( 5 ) 



sin 2 fl 



^-(w) = PU-M) - i&T^f • (6) 



These expressions are simple but approximate and hold at the Z pole, neglecting 
electromagnetic contributions, masses and terms proportional to g l v ; the full correct 
expressions can be found in Ref. [f|. 

Notice that, inserting the values of the coupling constants 

u,c,t = 1 _ 1 gin 2 Q u c,t = I / 7 n 

J y 2 3 2 

1 2 1 

„d,s,b _i_ c^n 2 /) „d,s,b 

9 = — sin v... q ' = 

jy 2 3 2 



one has 



sin 2 

P+ _._ + (uU,CC,ti) ~ ~ Q - 36 l + cos 2g ( 8 ) 



p + _._ + {dd 1 ss ) bb) ~ —0.17 



sin 2 1 



1 + cos 2 



Eq. (H) clearly shows the 6* dependence of p + _._ + (qq). In case of pure electro- 
magnetic interactions (-^/s -C M z ) one has exactly: 

Notice that Eqs. @ and @ have the same angular dependence, but a different sign 
for the coefficient in front, which is negative for the Z contribution. 

By using the above equations for p(qq) into Eq. ([I]) one obtains the most general 
expression of p(h) in terms of the qq spin state and the unknown fragmentation 
amplitudes. 
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Despite the ignorance of the fragmentation process some predictions can be made 
by considering the production of hadrons almost collinear with the parent jet: 
the qq —>■ h + X fragmentation is then essentially a cm. forward process and the 
unknown D amplitudes must satisfy the angular momentum conservation relation 

1 

where Oh is the angle between the hadron momentum, h = zq + p , and the quark 
momentum q, that is 

sin0 h ~%. (11) 

The bilinear combinations of fragmentation amplitudes contributing to p(h) are then 
not suppressed by powers of (p T / Z\fs) only if the exponent in Eq. fllPl) is zero, which 
greatly reduces the number of relevant helicity configurations. 

The fragmentation process is a parity conserving one and the fragmentation 
amplitudes must then also satisfy the forward parity relationship 

D -x h -x x ,- + = (-l) 5 * +5 * +A *- A * D KXx;+ _ . (12) 

Before presenting analytical and numerical results for the coherent quark frag- 
mentation let us remember that in case of incoherent single quark fragmentation 
Eq. ([[]) becomes 



ly h q,X,X x ,X q ,\' q 



where p(q) is the quark q helicity density matrix related to p(qq) by 

Pa v =£Pa„a -;A' (14) 



?' q> «' i 



In such a case angular momentum conservation for the collinear quark fragmen- 
tation requires X q = X h + X x ; the Standard Model computation of p(q) gives only 

non zero diagonal terms [p ++ (q) = p + _. + _(qq) and p (q) = p_ + ._ + (qq)], and one 

ends up with the usual probabilistic expression 

where D^\ h is the polarized fragmentation function of a q with helicity X q into a 
hadron h with helicity X h . Off-diagonal elements of p(h) are all zero. 
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2. e~e+ -> BX (5 fl = 1/2, p T /V^ -»• 0) 

Let us consider first the case in which h is a spin 1/2 baryon. It was shown in 
Ref. ||] that in such a case the coherent quark fragmentation only induces small 
corrections to the usual incoherent description 



P 



(B) 
P + -(B) 



N, 



O 



P+-,+-(qq) D* 



P- 



P 3 



Z\ s 



(16) 
(17) 



That is, the diagonal elements of p(B) are the same as those given by the usual 
probabilistic formula ([15]), with small corrections of the order of (p T / Zyfs) 2 , while 
off-diagonal elements are of the order {jp T /zyfs) and vanish in the p T /y/s — > limit. 

The matrix elements of p(B) are related to the longitudinal (P z ) and transverse 
(P y ) polarization of the baryon: 



2p 



++ 



1. 



-2 Imp H 



(18) 



Some data are available on A polarization, both longitudinal and transverse, from 
ALEPH Collaboration || and they do agree with the above equations. In particular 
the transverse polarization, at yfs = M z , p T ~ 0.5 GeV/c and z ~ 0.5 is indeed of 
the order 1%, as expected from Eq. (ITS 



3. e-e+ -> VX (S v = 1, p T /V* -»• 0) 

In case of final spin 1 vector mesons one has, always in the limit of small p T Jl[], 



Poo(V) 
Pn(V) 
Pi,-i(V) 



ly V q 



N, 



•(qq)Dq& + p-+-,-+{qq) D q~ 



V q 



(19) 

(20) 
(21) 



V q,X 



Again, the diagonal elements have the usual probabilistic expression; however, 
there is now an off-diagonal element, p x _ 1( which may survive even in the p T j \fs — > 
limit. In the sequel we shall concentrate on it. Let us first notice that, in the collinear 
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limit, one has 



d7 + = ElJW-l a = ^ ( 22 ) 

D V q 'X = Y,\D 10;+ -\ 2 = (23) 

= J2\Dn;- + \ 2 = D^\ (24) 

x 

with Z}^ = Z}^+ + + -D^+ 1 and N v = J2 q Dj. We also notice that the two frag- 
mentation amplitudes appearing in Eq. (21) are related by parity and their product 
is always real. p 00 and p 1 _ x can be measured through the angular distribution of 
two body decays of V. 

In order to give numerical estimates of p l _ l we make some plausible assumptions 

D h q A - = D^ = Q (25) 
D h q ;° + = a v q D h q X- (26) 

The first of these assumptions simply means that quarks with helicity 1/2 (—1/2) 
cannot fragment into vector mesons with helicity —1 (+1). This is true for valence 
quarks assuming vector meson wave functions with no orbital angular momentum, 
like in SU(6). The second assumption is also true in SU(Q) with = 1/2 for any 
valence q and V. Rather than taking = 1/2 we prefer to relate the value of 
to the value of p 00 (V) which can be or has been measured. In fact, always in the 
p T — > limit, one has || 

If is the same for all valence quarks in V (aY = a v ) one has, for the valence 
quark contribution: 

v _ Poo( V ) 



v = ruuvw . (2g) 
1-PooOO 1 ' 



Finally, one obtains 



P..-AV) - [1 - Po ,o(n gl^tff±M . (29) 

We shall now consider some specific cases in which we expect Eq. (p9|) to hold; 
let us remind once more that our conclusions apply to spin 1 vector mesons produced 
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in e~e + — > gg — ► V + X processes in the limit of small p T and large z, i.e., to vector 
mesons produced in two jet events (e~e + — > gg) and collinear with one of them 
(p T = 0), which is the jet generated by a quark which is a valence quark for the 
observed vector meson (large z). These conditions should be met more easily in the 
production of heavy vector mesons. 
One obtains 0: 

Pw (iT) ~ [1 - p 0i0 (ir)] p+_ ; _ + (66) (30) 
Pi,-i(D*) * [1 - Pofi (D*)} P+-,-+(cc) (31) 
p h _M ~ [l-p OiQ (0)]p + _ ; _ + (ss) (32) 

Pi-i(p) - g f 1 ~ ^o,o(p)] \P+--+(uu) + p+--+(dd)} (33) 

**) * \ [1 - Pb.oC^)] [>+-;-+(«*) + (34) 
Pi,-i(^°) ^ i[l-p 0i0 (^ )][p + - ; - + (rfrf) + p + - ; - + ( S s)]. (35) 

Eqs. (BCT)-(B5p show how the value of p 1 _ 1 (V) are simply related to the off- 
diagonal helicity density matrix element p+- ; -+(gg) of the gg pair created in the 
elementary e~e + — > qq process; such off-diagonal elements would not appear in the 
incoherent independent fragmentation of a single quark, yielding p 1 _i(V) = 0. 

By inserting into the above equations the value of p 00 when available || and the 
expressions of p + „._ + , Eq. (8), one has: 

/ „n / s sin 2 9 , , 

Pi-i(B ) ~ -(0.109 ±0.015) i + cos2& (36) 

Pi-i(D*) ~ -(0.216 ±0.007) 1 ^ 2g (37) 

/ n \ sin 2 . , 

p w (0) * -(0.078 ± 0.014) - (38) 

Pi,-i(K*°) * -0-170 [l-p , (iT°)] 1 ^ 2g - (39) 

Finally, in case one collects all meson produced at different angles in the full available 
range (say a < < tt — a, |cos#| < cos a) an average should be taken in 0, 
weighting the different values of p x _i(9) with the cross-section for the e~e + — > V+X 
process; this gives ||: 

3 — cos 2 CY 

(Pi,-i(£*)W- a] * -(0.109 ± 0.015) — — (40) 

6 + cos z a 
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•(0.216 ±0.007) 



•(0.078 ±0.014) 



(Pl,-l(0)) [«,*-«] ^ 

These results have to be compared with data 



cos 2 a 



3 ± cos 2 a 
3 — cos 2 a 



3 ± cos 2 a 
0.170 [l-p JK*)] 6 



cos 2 a 



3 ± cos 2 a 



(41) 
(42) 
(43) 



Pi,-i(D*) = 


-0.039 ±0.016 


for 


z > 0.5 


cos a = 


0.9 


(44) 


Pi-x{4>) = 


-0.110 ±0.070 


for 


z > 0.7 


cos a = 


0.9 


(45) 


Pi,-i(K*°) = 


-0.090 ±0.030 


for 


z > 0.3 


cos a = 


0.9 


(46) 



which shows a good qualitative agreement with the theoretical predictions. We 
notice that while the mere fact that pi~i differs from zero is due to a coherent 
fragmentation of the qq pair, the actual numerical values depend on the Standard 
Model coupling constants; for example, pi-i would be positive at smaller energies, 
at which the one gamma exchange dominates, while it is negative at LEP energy 
where the one Z exchange dominates. has also a peculiar dependence on the 

meson production angle, being small at small and large angles and maximum at 
9 = 7r/2. Such angular dependence has been tested in case of K*° production and 
indeed one has 0, in agreement with Eqs. ( |35"D and (§), 



' Pl,-1 ' 






. 1 ~Poo. 


|cos0|<O.5 


. 1 ~Poo. 



1.5 ±0.7 



(47) 



J |cos6»|>0.5 



4. Diagonal elements of p(V) and P v = V/(V + P) 

Let us consider now the diagonal element poo(^) ~ for which the probabilistic in- 
terpretation, Eq. ( |I9|) holds - together with the production of pseudoscalar mesons; 
that is, we consider the production of the pseudoscalar mesons P = K,D,B and 
the corresponding vector mesons V = K*, D*, B*\ data are available on poo(^) and 
the ratio of vector to vector + pseudoscalar mesons, P v =V/(V + P) || 

We denote by Pg the probability that the fragmenting quark produces a meson 
with spin S and helicity A and consider only the production of vector and pseu- 
dovector mesons. Then we have || |9|: 

Poo(V) = -J± o0 P y = P 1 ±1 + P 1 ° (48) 
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with P^ 1 = Pi + Pf 1 and P^ 1 + Ff + P ° = 1. 

In terms of the fragmentation functions this reads: 

P p P y '° P^ 1 + D^- 1 

pO = pO _ p±l = Z * ( AO) 

D v + D P 1 D v + D P 1 D v + D P y ' 

q 1 q q 1 q q 1 q 

Notice that, by parity invariance, the above quantities are independent of the quark 
helicity. 

Statistical spin counting would give 

p±! = 0.5 Pf = 0.25 P ° = 0.25 (50) 

that is 

Poo = \ Pv = l (51) 
so that the vector meson alignment is zero: 

A = \ (3p 00 - 1) = (52) 
From data [0 one obtains [0 

Ff\K*) = 0.34 ±0.06 P?(K*) = 0.41 ± 0.07 Pq{K) = 0.25 ± 0.10 (53) 

Pf^CD*) = 0.34 ± 0.04 Pi(D*) = 0.23 ± 0.03 P °(P) = 0.43 ± 0.06 (54) 
P±\B*) = 0.49 ± 0.09 Pi(B*) = 0.27 ± 0.08 P °(P) = 0.24 ± 0.09 (55) 

The simultaneous measurements of the ratio of vector to vector + pseudovector 
mesons and Poo(^) supply basic information on the fragmentation of quarks which 
does not depend on the helicity of the quark, but on the spin and helicity of the 
final meson. The data available for K, D and B mesons show clear deviations from 
simple statistical spin counting; such information could be of crucial importance for 
the correct formulation of quark fragmentation Monte Carlo programs, which at the 
moment widely assume simple relative statistical probabilities. 

Let us consider our results, Eqs. fl53"|)-(|55l). For strange mesons the data agree 
with spin counting in the amount of K versus K*, but, among vector mesons, helicity 
zero states seem to be favoured; these are in absolute the most abundantly produced, 
P®(K*) = 0.41. For charmed mesons results differ from the spin counting values 
(|50|), suggesting a prevalence of pseudoscalar states, Pq(D) = 0.43, and, among 
vector mesons, of helicity states. The heavy 6-mesons, instead, are produced in 
good agreement with statistical spin counting rules, as one expects. 
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5. pi,_i(V) in other processes and conclusions 

The results discussed here are encouraging; indeed measurements of off-diagonal 
and diagonal elements of p(V) give valuable information on the hadronization pro- 
cess and test the underlying elementary dynamics. It would be very helpful to have 
more and more detailed data, possibly with a selection of final hadrons with the 
required features for our results to hold. 

It would be interesting to test the coherent fragmentation of quarks in other 
processes ||1Q|| , like 77 — > VX, pp — > D*X and 7p — > VX or 7*p — > VX. The 
first two processes are similar to e~e + — > VX in that a qq pair is created which 
then fragments coherently into the observed vector meson; one assumes that the 
dominating elementary process in pp — > D*X is gg — > cc. In both these cases one 
has for p + _. _+(<?<?) the same value as in Eq. so that one expects a positive value 
<>f>. lGO- 

In the case of the real photo-production of vector mesons the quark fragmentation 
is in general a more complicated interaction of the struck quark with the remnants 
of the proton and it might be more difficult to obtain numerical predictions. How- 
ever, if one observes D* mesons one can assume or select kinematical regions for 
which the underlying elementary interaction is 7^ — > cc: again, one would have the 
same p + _._ + (cc) as in Eq. (|9]), and one would expect a positive value of p 1 _ l {D*). 
Similarly for the production of or B* . 

The production of vector mesons like D* or B* in DIS is even more interesting; 
the polarization of the virtual photon depends on x and Q 2 and so does the value of 
. An eventual dependence of p l _i(D*) on x and Q 2 would then be an 
unambigous test of the hadronization mechanism and the elementary interaction. A 
similar situation can be obtained by considering e + e~ or 77 processes with polarized 
initial particles: in such cases the value of p + _. _+(<?<?) strongly depends on the initial 
spin states and should change the measured value of p 1 |p"T] . 

Similar considerations hold for the diagonal elements of p(V) and for P v ; the 
degree of universality of quark fragmentation could be tested by studying these 
same quantities in other processes, like the ones mentioned above. It would also be 
interesting to compare data on the production of spin 1/2 and spin 3/2 baryons. 

To conclude, some non perturbative aspects of strong interactions can only be 
tackled by gathering experimental information and looking for patterns and regular- 
ities which might allow the formulation of correct phenomenological models. Spin 
dependent quantities are still in a first stage of consideration and development, so 
that even qualitative studies are meaningful; for example, it would indeed be in- 
teresting to perform the simple tests of coherent fragmentation effects suggested 
here. 



P+- 



cc 



m 
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